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SONAR  BEAMFORMING  -  AN  OVERVIEW  OF  ITS  HISTORY  AND 

STATUS 


INTRODUCTION 


This  report  introduces  the  subject  of  sonar  beamforming.  It  is  assumed 
that  the  reader  has  little  prior  knowledge  of  the  topic.  A  narrative  approach  has 
been  taken  that  follows  a  historical  thread  from  the  early  days  of  beamforming  to 
the  present.  Numerous  references  are  given  for  papers  and  texts  should  the 
reader  desire  to  further  explore  the  subject.  Only  unclassified  material  is 
presented  under  this  cover;  related  classified  topics  will  be  treated  in  a 
confidential  addendum. 


DIRECTIONAL  ANTENNAS 


WHAT  ARE  THEY? 

In  general,  an  antenna  is  a  transducer  (or  group  of  transducers)  between 
an  electrical  waveform  in  a  circuit  and  a  field  of  energy.  For  example,  an 
antenna  might  sense  electromagnetic  energy  (radio  waves)  or  mechanical 
energy  (sound  pressure  fluctuations).  Suppose  that  an  antenna  is  placed  at  the 
center  of  a  sphere  of  infinite  radius  and  that  we  observe  its  output  as  a  radiating 
source  of  energy  is  moved  along  the  sphere's  surface.  If  the  output  power  of  the 
antenna  remains  invariant  to  the  position  (arrival  direction)  of  the  source,  we 
say  that  the  antenna  is  omnidirectional.  If,  however,  the  output  power  varies 
with  the  arrival  direction  of  the  source,  then  we  say  that  the  antenna  is 
directional.  An  ideal  directional  antenna  would  have  zero  output  power  (i.e., 
infinite  attenuation)  when  the  source  was  arriving  from  all  directions  but  one;  the 
direction  to  the  source  when  the  antenna  output  is  maximum  is  called  the  steer 
(or  look)  direction,  in  practice,  actual  directional  antennas  maintain  nearly 
constant  response  to  signal  fields  arriving  within  a  finite  angular  sector  atout 
the  steer  direction  and  provide  only  finite  amounts  of  attenuation  to  those 
arriving  outside  that  sector. 


WHY  DO  WE  NEED  THEM? 

There  are  two  reasons  why  we  use  directional  antennas.  First,  we  may 
not  always  know  a  priori  the  direction  of  the  source  we  are  trying  to  receive.  An 
omnidirectional  antenna  might  serve  to  detect  the  presence  of  the  source  in 
such  a  situation,  but  would  give  no  indication  of  where  it  is  located.  A 
directional  antenna,  however,  allows  determination  of  the  source  arrival  angle 
because  its  output  power  is  maximum  when  steered  at  the  source. 
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The  second  reason  for  using  a  directional  antenna  is  to  attenuate  signals 
that  arrive  from  directions  other  than  the  steer  direction.  Such  signals  interfere 
with  our  ability  to  observe  a  waveform  arriving  from  the  steer  direction  and 
hence  act  as  noise.  By  attenuating  the  contribution  from  interfering  sources,  a 
directional  antenna  improves  the  signaMo-noise  ratio  (SNR),  which  is  the  ratio 
of  signal  power  to  noise  power. 


CONTINUOUS  APERTURES  VERSUS  ARRAYS 

An  antenna's  aperture  is  the  region  in  space  within  which  it  senses 
energy.  Given  an  aperture,  we  can  design  a  directional  antenna  in  one  of  two 
ways.  We  can  either  populate  the  aperture  with  a  continuous  sensor  or  we  can 
sample  the  aperture  with  a  finite  number  of  sensors.  An  example  of  a 
continuous  aperture  is  the  parabolic  dish  shown  in  figure  1 ;  this  particular  dish 
is  used  to  make  radio  astronomy  observations.  As  with  ail  continuous  aperture 


Rgure  1 .  Example  of  Continuous  Aperture 


antennas,  the  steer  direction  of  the  dish  is  changed  by  mechanically  changing 
the  orientation  of  the  device.  Although  it  may  pose  no  serious  problem  to 
mechanically  steer  an  antenna  for  some  applications,  such  as  radio  astronomy, 
other  applications  may  preclude  mechanical  steering.  For  example,  scanning 
the  steer  direction  of  a  large,  shipboard  air  surveillance  antenna  may  be 
required  at  rates  too  fast  to  achieve  mechanically.  In  this  case,  we  employ  the 
sampled  aperture  approach.  The  set  of  sensors  used  to  sample  the  aperture  is 
called  an  array,  and  as  we  shall  see  in  a  following  section,  we  replace 
mechanical  steering  with  an  electrical  steering  process  called  beamtorming. 
Figure  2  shows  a  full-scale  test  fixture  for  two  arrays  that  will  ultimately  be 


installed  in  the  bow  area  of  a  submarine. 


Figure  2.  Examples  of  Bow-Mounted  Submarine  Arrays 


BEAM  PATTERNS 

in  order  to  graphically  represent  the  directional  properties  of  an  antenna, 
we  use  a  plot  called  a  beam  pattern.  We  obtain  such  a  pattern  by  plotting 
output  power  as  a  function  of  the  arrival  angle  of  a  plane  wave  source.  Thus,  a 
beam  pattern  is  a  plot  in  three  dimensions.  A  sample  pattern  is  shown  in  figure 
3. 


We  see  that  the  pattern  has  one  lobe  that  is  taller  than  ail  the  others.  This 
is  called  the  main  lobe,  and  its  angular  extent  is  called  the  beamwidth.  The  axis 
through  the  center  of  the  main  lobe  is  called  the  maximum  response  axis  (MRA). 
The  angular  sector  outside  of  the  main  lobe,  called  the  sidelobe  region, 
contains  all  the  lower  lobes  of  the  pattern,  it  is  this  region  which  provides 
attenuation  of  signals  arriving  from  directions  other  than  the  MRA.  It  is  standard 
practice  to  plot  beam  patterns  in  units  of  decibels  (dB) .  If  6  and  ^  are  polar 
angles  and  5  (e,i|>)  is  the  output  power  when  the  source  arrives  from  direction 
(e,<|>),  then  the  decibel  value  of  the  beam  pattern  is  given  by 

B(e,(|i)  =  ioiog,„[b(e.«] 

where 

b(e,(|))  =  B(0.<|i)/6(9,,(|i,)  (2, 
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Figure  3.  Sample  Beam  Pattern 


is  a  scaled  version  of  b  (0,<|>),  which  has  been  normalized  to  a  maximum  value 
of  unity,  and  where  (0s.<|>s)  is  the  steer  direction  of  the  MRA.  It  is  also  standard 
practice  to  define  beamwidth  as  the  number  of  degrees  between  the  3  dB  down 
points  of  the  mainlobe. 


ARRAY  GAIN 

One  might  wonder  to  what  extent  a  directional  antenna  will  improve  SNR 
over  an  omnidirectional  antenna.  To  characterize  this  improvement,  we  define 
a  quantity  called  array  gain  (AG)  as  the  ratio  of  the  output  SNR  of  a  directional 
antenna  to  the  output  SNR  of  an  omnidirectional  antenna  in  the  same  signal- 
plus-noise  field;  i.e.. 


AG  = 


directional 

(SNR)  , 

'  '  omnidiractional 


(3) 


As  with  beam  patterns,  it  is  common  to  express  array  gain  in  decibels  (dB)  by 
computingtO  times  logio(AG). 


DIRECTIVITY  INDEX 


The  array  gain  of  a  directional  antenna  will  vary  with  the  directional 
nature  of  the  noise  field  ,  as  well  as  with  the  coherence  of  the  signal  field  across 
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the  aperture.  For  better  or  worse  (probably  the  latter),  many  designers  and 
users  of  antennas  prefer  to  assign  a  single  number  :o  an  antenna  that 
characterizes  its  directional  behavior.  The  number  is  called  directivity  index 
(Dl)  and  is  defined  as  the  array  gain  of  the  antenna  for  the  special  case  in 
which  (1)  a  perfectly  coherent  plane  wave  signal  arrives  on  the  MRA,  and  (2) 
the  noise  field  is  isotropic  (i.e.,  the  same  amount  of  noise  arrives  from  all 
directions).  It  can  be  showni  that  when  the  response  of  a  directional  antenna  to 
a  plane  wave  signal  arriving  on  it's  MRA  equals  the  response  of  an 
omnidirectional  antenna  to  the  same  signal,  the  Dl  is  given  by 

^  Volume  of  Unit  Sphere 
Volume  of  b(8,<j)) 

Recall  that  b(e.<t>)  is  the  beam  pattern  whose  maximum  value  on  the  MRA  has 
been  normalized  to  unity.  Thus,  the  volume  of  b(6,()>)  ^  volume  of  unit  sphere, 
and  so.  from  equation  (4).  we  have  Dl  ^  1.  Also  from  equation  (4).  we  see  that 
as  the  volume  contained  within  a  beam  pattern  is  decreased  (e.g..  by  narrowing 
the  main  lobe  and/or  lowering  the  side  lobes),  Dl  is  increased. 


SINGLE  BEAM  SYSTEMS  (1940's  •  1960's) 


MECHANICAL  STEERING  OF  CONTINUOUS  APERTURES 

The  only  way  to  change  the  look  direction  of  a  continuous  aperture 
antenna  is  to  mechanically  move  its  MRA.  During  World  War  II.  U.  S. 
submarines  were  equipped  with  the  JP  sonar  that  used  a  continuous  aperture 
in  the  form  of  a  rotatable  horizontal  line  hydrophone.  As  shown  in  figure  4,  the 
operator  steered  the  line  hydrophone  via  a  servo  system  to  bearings  of  interest. 
A  major  iimitation  of  such  mechanically  steered  antennas  is  size.  It  is  not  easy 
to  fabricate  large  hydrophones,  and  even  if  it  were,  the  size  could  never  exceed 
the  ship's  width.  To  obtain  larger  apertures  it  became  necessary  to  construct 
arrays  of  many  small  hydrophones  and  then  to  steer  them  electrically  as 
described  in  the  next  section. 


BEAMFORMING  (ELECTRICAL  STEERING) 

The  principle  behind  beamforming  is  really  quite  simple.  Consider  a 
plane  wave  signal  field  impinging  upon  an  array  of  hydrophones.  An 
examination  of  the  output  waveforms  from  the  hydrophones  would  show  that 
they  are  delayed  versions  of  each  other,  with  the  amount  of  delay  between  any 
given  pair  equal  to  the  time  required  for  the  signal  to  propagate  between  the 
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Figure  4.  Mechanically  Steered  JP  Sonar 


pair.  For  example,  consider  a  signal  field  arriving  at  a  line  array  of  M 
equispaced  hydrophones  as  shown  in  figure  5.  For  the  geometry  shown,  the 


Figure  5.  Geometry  For  Signal  Reid  Arriving  At  Line  Array 


signal  field  first  arrives  at  hydrophone  1 .  then  at  2. 3,  . .  M.  Thus,  the  output 
waveforms  from  the  hydrophones  would  appear  as  illustrated  in  figure  6. 
in  addition  to  the  signal,  each  hydrophone  output  will  also  contain  undesired 
noise.  Because  the  noise  tends  to  be  independent  from  phone  to  phone,  time¬ 
aligning  the  signal  waveforms  and  then  summing  will  cause  the  signals  to  add 
coherently  and  the  noises  incoherently,  thus  improving  SNR.  This  process  of 
time-alignment  followed  by  summation  is  called  conventional  time  domain 
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n  •  hydreplWM  Indsx 


Figure  6.  Output  Waveforms  From  Hydrophones  Due  to  Signal  Field 
beamforming  and  is  shown  in  block  diagram  form  in  figure  7. 


Figure  7.  Conventional  Time  Domain  Beamformer  Block  Diagram 


in  some  applications,  hydrophone  outputs  are  passed  through 
narrowband  filters  before  beamforming.  This  filtering  may  be  done  because  it 
is  known  a  priori  that  the  signal  is  narrowband,  or  perhaps  because 
implementation  is  preferred  in  the  frequency  domain  via  a  fast  Fourier  transform 
(FFT).  Regardless  of  the  motivation,  once  the  hydrophone  outputs  have  been 
narrowband  filtered,  it  is  possible  to  approximate  the  time  delays  needed  for 
beamformation  by  phase  shifts.  If  y(fo.t)  is  the  complex  envelope  out  of  a 
narrowband  filter  having  bandwidth  B  and  center  frequency  fo.  then  it  can  be 
shown  that  a  good  approximation  to  a  delayed  version  of  y(fo,t)  is  given  by 

y(fo*t -T)  5exp(-j27cfoT)  y(fo,t) ,  BT<0. 1  (5) 


A  block  diagram  of  a  conventional  frequency  domain  beamformer  is  shown  in 
figure  8. 
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Rgure  8.  Conventional  Frequency  Domain  Beamformer  Block  Diagram 


From  the  above  discussion,  we  see  that  changing  the  steer  direction  of  a 
beamformer  requires  changing  the  time  delays  (or  phase  shifts)  applied  to  the 
hydrophone  outputs.  It  was  the  Germans  who  found  a  practical  way  to 
accomplish  this  during  World  War  II.  They  developed  the  so-called 
compensator  system^  shown  in  figure  9.  The  system  consists  of  two  plates  - 


Figure  9.  Compensator  Diagram 


one  is  fixed  and  the  other  rotates.  The  fixed  plate  is  covered  with  straight, 
narrow,  evenly  spaced  strips  of  metal,  each  of  which  is  wired  to  the  input  of  a 
tapped  delay  line.  The  movable  plate  has  brushes  that  slide  over  the  fixed 
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plate.  The  brushes  are  arranged  in  the  same  geometric  shape  as  the  array,  and 
each  brush  is  wired  to  its  corresponding  hydrophone.  The  operator  is  given  a 
handwheel  that  controls  the  rotation  of  the  movable  plate.  As  the  operator 
steers  the  handwheel,  the  brushes  on  the  movable  plate  are  connected  to  the 
proper  taps  on  the  delay  line,  and  the  superposition  (sum)  of  the  delayed 
outputs  appears  at  the  output  of  the  delay  line.  Manual  operation  of  a 
compensator-based  system  is  illustrated  in  figure  10. 


Figure  10.  Manual  Operation  of  a  Compensator-Based  System 


The  advent  of  beamforming  by  means  of  compensators  gave  sonar 
operators  the  capability  to  manually  scan  a  single  beam  in  any  direction  without 
requiring  movement  of  the  array  itself.  Next,  as  shown  in  figure  1 1 ,  it  was  a 
simple  step  to  build  systems  with  two  compensators,  one  for  manual  use  by  the 
operator  and  the  other  connected  to  a  motor  that  repeatedly  scanned  its  beam 
through  all  bearings.  The  detected  output  of  the  electrically  scanned  beam  was 
used  to  intensity  modulate  a  strip  chart  recorder  in  a  bearing  versus  time  format, 
thus  providing  the  operator  with  a  time  history  of  traces  from  all  detected 
contacts.  The  operator  could  then  use  the  nianually  scanned  beam  to 
investigate  contacts  of  interest  for  classification  purposes. 

Electrically  scanned  sonars  were  a  major  improvement  over 
mechanically  scanned,  continuous  aperture  antennas.  However,  they  still 
suffered  from  a  serious  drawback  called  scanning  loss.  A  scanned  beam 
spends  only  a  small  fraction  of  its  time  pointed  at  a  given  target.  In  fact,  if  the 
beamwidth  is  a  degrees,  then  the  beam  is  trained  on  a  contact  only  a/360 
percent  of  the  timel  To  eliminate  scanning  loss,  the  sonar  community  turned  its 
attention  to  the  development  of  multibeam  systems,  which  is  the  subject  of  the 
next  section. 
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Figure  1 1 .  Compensator-Based  System  With  Continuous  and  Manual 

Scanning 


MULTIBEAM  SYSTEMS  (1960's  -  PRESENT) 


DIGITAL  MULTIBEAM  STEERING  (DIMUS) 

To  eliminate  scanning  loss,  beamformers  were  developed  that 
simultaneously  form  multiple  beams  pointed  to  fixed  directions  in  space.  These 
beams  collectively  provide  continuous  coverage  in  time  over  any  desired 
spatial  sector.  Figure  12  illustrates  a  case  in  which  the  multiple  beams  of  an 
array  are  providing  complete  coverage  in  azimuth,  and  limited  coverage  in 
vertical  angle  (to  reduce  clutter  in  the  figure,  a  complete  set  of  vertical  beams  is 
shown  at  only  one  azimuth).  This  approach  ensures  that  one  beam  is  always 
pointed  at  any  given  target  in  the  sector  being  covered. 

One  of  the  first  multibeam  systems  was  built  by  Dr.  Victor  Anderson  of  the 
Marine  Physical  Laboratory,  San  Diego,  California,  in  the  late  ISSO's^.  Digital 
circuits  were  also  becoming  available  at  that  time,  and  Anderson  took 
advantage  of  this  new  technology  to  implement  the  delay  function  required  for 
beamforming.  As  shown  in  figure  13,  he  converted  the  outputs  of  hydrophones 
into  bilevel  waveforms  by  passing  the  signsds  through  clipper  amplifiers,  and 
then  he  fed  the  bilevel  waveforms  into  shift  registers  which  senred  as  tapped 
delay  lines.  Outputs  from  the  appropriate  shift  register  stages  were  fed  into 
analog  summing  amplifiers  to  form  the  multiple  beams.  Eventually,  the  analog 
summing  amplifiers  were  replaced  by  digital  accumulators,  yielding  an  all 
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■■V 

Figure  12.  Illustration  of  Multiple  Beams 


digital  beamformer  with  no  moving  parts  and  no  scanning  loss.  This  was  a 
major  advance  in  beamformer  technology. 


Rgure  1 3.  Schematic  Diagram  Of  The  DIMUS  Technique 


LINEAR  MULTIBEAM  SYSTEMS 

While  the  DIMUS  beamformer  was  a  great  step  forward,  it  still  had  one 
drawback.  Clipping  the  hydrophone  outputs  prior  to  beamforming  introduced 
spectral  distortion  that  was  not  acceptable  for  some  applications.  Fortunately, 
the  rapid  development  of  digital  integrated  circuits  during  the  1970's  and  80’s 
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made  it  possible  to  replace  the  clipper  amplifiers  with  analog-to-digital 
converters,  and  the  one  bit  shift  registers  with  multibit  memories.  Thus,  it 
became  possible  to  maintain  the  linearity  of  hydrophone  waveforms  throughout 
the  time  domain,  beamforming  process,  in  a  similar  fashion,  the  advent  of 
digital  hardware  to  compute  FFTs  allowed  formation  of  multiple,  linear  beams 
in  the  frequency  domain. 


ADAPTIVE  BEAMFORMING  (ABF) 


STATIC  AMPLITUDE  SHADING 

Before  discussing  adaptive  beamforming,  we  shall  consider  a  related 
topic  called  shading  of  an  array.  Early  array  designers  looked  at  beam  patterns 
and  wondered  if  it  might  be  possible  to  modify  the  beamforming  process  to  yield 
lower  side  lobe  levels.  The  answer  is  yes.  if  a  widening  of  the  pattern's  main 
lobe  is  acceptable.  The  technique,  as  shown  in  figure  14.  involves  multiplying 
each  hydrophone  output  by  a  properly  chosen  real  gain  prior  to  either 
conventional  time  or  frequency  domain  beamformation. 


Xi(t)  »  Ai  1  '  » 

Conventional 

Time 

Domain  or 

- ►z(t) 

Frequency 

Domain 

Beamfonner 

Xj^(t)  ^  Am  1  » 

Rgure  14.  Conventional  Beamformer  With  Amplitude  Shading 


The  set  of  gains,  {Aj}.  are  called  amplitude  shading  weights.  An  array  in  which 
all  shading  weights  are  equal  is  said  to  be  uniformly  shaded  and  is  equivalent 
to  having  no  shading  at  all.  Various  criteria  have  been  used  to  "properly 
choose"  a  set  of  nonuniform  shading  weights.  For  example,  one  might  ask, 
"What  set  of  weights  will  give  the  narrowest  main  lobe  while  keeping  all  side 
lobe  peaks  at  a  specified  level?"  A  procedure  for  computing  such  a  set  of 
weights  was  derived  by  Dolph  using  Chebyshev  polynomials,  resulting  in  what 
is  known  as  Dolph-Chebyshev  shading.  Then  again,  one  might  not  desire  all 
side  lobe  peaks  to  be  at  the  same  level,  but  rather  to  drop  off  at  a  prescribed 
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rate.  This  can  be  achieved  using  cos°(x)  shading,  where  a  controls  the  rate  of 
fall-off.  These  are  only  two  examples  from  the  wide  variety  of  shading  schemes 
that  have  been  developed  over  the  years. 


ABF  AS  A  GENERALIZATION  OF  SHADING 

Three  points  discussed  in  the  preceding  section  are  particularly 
noteworthy.  First,  criteria  for  determination  of  the  shading  weights  were  based 
solely  upon  the  array's  response  to  the  signal  field;  no  mention  was  made  of  the 
response  to  the  noise  field.  Second,  the  shading  weights  were  constrained  to 
be  real,  even  for  frequency  domain  beamforming.  And  third,  no  mention  was 
ever  made  of  changing  the  shading  weights  once  they  were  computed.  In  the 
early  1960's  researchers  took  a  more  general  approach  to  array  processing. 
They  included  signal  aojl  noise  properties  in  the  derivation  of  shading  weights, 
and  because  these  properties  could  change  with  time,  they  decided  that  the 
shading  weights  should  be  adapted  accordingly  (hence  the  term  adaptive 
beamforming ).  Although  any  discussion  of  ABF  can  be  conducted  in  either  the 
time  or  frequency  domain,  we  shall  use  with  the  latter  as  has  been  done  in  most 
of  the  literature  and  in  most  actual  implementations.  Because  the  inputs  to  a 
frequency  domain  beamformer  may  be  represented  by  complex  numbers  (e.g., 
FFT  coefficients  or  complex  envelopes),  a  further  generalization  will  be  made  by 
allowing  the  shading  weights  to  be  complex  also.  Thus,  our  ABF  block  diagram 
will  be  as  shown  in  figure  15,  where  the  {Xj}  and  {Wj}  are  complex  and  the  {WJ 
are  adjustable. 


Figure  15.  Frequency  Domain  ABF  Block  Diagram 


OPTIMALITY  CRITERIA  FOR  ABF 


By  the  early  1960's.  single  channel,  optimum  filter  theory  for  signal 
detection  and  estimation  was  fairly  well  established.  Van  VIeck  and  Middleton^ 
had  derived  the  "matched  filter"  that  maximized  SNR  for  detection,  and  Wiener^ 
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had  derived  the  filter  for  estimation  that  minimized  mean  square  error  (MMSE). 
Subsequent  researchers  derived  filters  based  on  other  criteria  of  optimality, 
e.g.,  maximum  likelihood  ratio  for  detection  and  maximum  likelihood  function  for 
estimation.  It  was  only  natural  that  these  results  be  extended  to  the 
multichannel  application  of  beamforming.  Bryn^  did  so  for  the  maximum 
likelihood  ratio  case.  Mermoz^  for  maximum  SNR,  Burg^  for  MMSE.  and 
Darlington^  and  Levin^Q  for  minimum  variance  distortionless  response  (MVDR). 
Under  the  MVDR  criterion,  a  set  of  filter  (shading)  weights  is  found  that 
minimizes  power  (variance)  out  of  the  beamformer  while  constraining  the 
beamformer  to  pass  the  signal  with  fixed  gain  and  no  spectral  distortion. 

In  1966  Van  Trees^^  published  a  significant  paper  showing  that  ail  the 
various  optimum  array  processors  mentioned  above  could  be  represented  by 
the  form  shown  in  figure  16,  i.e.,  an  MVDR  beamformer  followed  by  a  single 
filter  whose  transfer  function  depends  upon  the  criterion  of  optimality. 


Figure  16.  Form  Of  All  Optimum  Array  Processors 


Because  the  MVDR  beamformer  is  common  to  all  optimum  array  processing 
schemes,  it  received  much  attention  from  researchers  during  the  1970's  and 
80's.  From  this  point  on,  whenever  we  refer  to  an  adaptive  beamformer,  we 
shall  assume  its  design  is  based  upon  the  MVDR  criterion  unless  otherwise 
stated. 


NULL  STEERING  INTERPRETATION  OF  ABF  PERFORMANCE 

As  noted  earlier,  computations  of  static  shading  weights  for  conventional 
beamformers  are  based  upon  criteria  that  are  independent  of  the  noise  field  in 
which  the  array  is  operating.  Thus,  it  is  possible  to  have  a  side  lobe  pointing 
toward  a  source  of  directional  interference.  This  situation  is  shown  in  figure  17. 
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Theta  (degrees) 

Figure  17.  Conventional  Beam  Pattern  With  Directional  Interference 


An  adaptive  beamformer,  on  the  other  hand,  will  adjust  its  weights  to  minimize 
output  power  while  maintaining  unit  gain  for  signals  arriving  from  the  steer 
direction.  The  way  it  minimizes  output  power  is  by  attenuating  its  response  to 
the  directional  interference.  In  terms  of  beam  patterns,  it  adjusts  its  side  lobe 
response  such  that  a  notch  or  null  is  positioned  at  the  arrival  angle  of  the 
interference  as  shown  in  figure  18. 


Theta  (degrees) 

Figure  18.  Adaptive  Beam  Pattern  With  Directional  Interference 


If  the  arrival  angle  of  the  interference  changes,  the  adaptive  beamformer 
automatically  "steers"  its  null  so  that  it  is  always  pointed  at  the  intenerence  - 
hence,  the  interpretation  of  an  adaptive  beamformer  as  a  null  steering  device. 
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ELEMENT  SPACE  VERSUS  BEAM  SPACE  ABF 


In  some  applications,  it  may  be  much  easier  to  access  the  outputs  of 
linear,  conventional,  preformed  beams  rather  than  the  outputs  of  individual 
array  elements.  In  such  a  case,  the  set  of  preformed  beams  corresponds  to 
outputs  from  an  array  of  directional  elements  (as  opposed  to  the 
omnidirectional  elements  of  the  original  array)  and  can  be  used  as  inputs  to  an 
MVDR  beamformer  as  shown  in  figure  19. 


Figure  19.  Element  Space  And  Beam  Space  Processors 


When  element  data  are  the  input  to  an  MVDR  beamformer,  we  use  the  term 
element  space  processor;  when  preformed  beam  data  are  the  input,  we  use  the 
term  beam  space  processor^^.  Gray^^  has  shown  that  as  long  as  the  number 
of  independent  beams  equals  the  number  of  array  elements,  the  beam  space 
processor  indeed  achieves  the  same  array  gain  as  the  element  space 
processor. 


ADAPTIVE  IMPLEMENTATIONS  OF  MVDR  BEAMFORMER 

Central  to  implementation  of  an  MVDR  beamformer  is  the  requirement  to 
solve  the  following  matrix  equation  for  the  optimum  weights: 

(6) 


where  w  is  the  optimum  weight  vector,  R,  is  the  cross-spectral  density  matrix  of 
the  input  data,  and  d  is  a  steering  vector  to  point  the  beam  in  some  desired  look 
direction.  There  are  basically  five  ways  in  which  this  equation  can  be 
implemented.  They  will  be  described  in  the  order  of  their  historical  evolution. 


Method  1 :  Estimate  R^^and  Invert 

This  could  be  called  the  brute  force  approach.  The  input  data  are  used 
either  recursively  or  in  a  sequential  batch  fashion  to  estimate  matrix  .  Each 
time  an  updated  estimate  of  is  generated,  it  is  directly  inverted  to  solve  for  w. 
Because  the  number  of  arithmetic  operations  required  to  invert  an  M  x  M  matrix 
is  0(M3),  where  0(  • )  denotes  "of  order,"  the  amount  of  computation  required 
can  be  substantial  for  large  M.  Furthermore,  numerical  instability  can  occur 
when  inverting  large  matrices  with  finite  precision  arithmetic.  Such  drawbacks 
are  the  reason  this  method  has  never  received  serious  consideration  for  real¬ 
time  implementation. 

Method  2:  Gradient  Descent 

Widrow^^  applied  the  method  of  gradient  descent  to  iteratively  solve  for 
w.  His  result,  known  as  the  least  mean  squares  (LMS)  algorithm,  never 
requires  the  estimation  of  R^  and  hence  never  requires  its  inversion.  The 
resulting  computational  requirement  is  0(M),  which  makes  it  extremely  attractive 
for  real-time  implementation.  It  is  particulartly  well  suited  for  running  on  a 
parallel  processor  of  either  the  systolic  array  type  or  the  single  instruction, 
multiple  data  (SIMD)  type.  Its  main  drawback  is  a  much  slower  convergence 
time  (the  number  of  iterations  required  to  get  sufficiently  "close"  to  the  optimum 
solution)  than  that  of  methods  3-5. 

Method  3:  Recursively  Estimate  •'* 

In  the  early  1970's  researchers  focused  on  ways  to  improve  speed  of 
convergence.  In  their  approach, based  upon  earlier  recursive  least  squares 
work  done  by  Plackett^®,  new  data  are  used  to  recursively  estimate  Rx rather 
than  Rx.  This  technique  saves  inverting  Rx,  and  the  method  indeed  converges 
faster  than  with  the  LMS  algorithm.  However,  the  price  paid  is  that  its 
computational  complexity  is  0(M2),  and,  unfortunately,  the  method  is  not  well 
suited  for  implementation  on  parallel  processors. 

Method  4:  Recursively  Estimate  Cholesky  Factor  and  Backsolve 

Because  Rx  is  a  Hermitian,  positive  definite  matrix,  it  can  always  be 
written  in  factored  form  as^^ 

R^-CC 

where  C  is  a  lower  triangular  matrix  known  as  the  Cholesky  factor.  Owsley^ ^ 
has  shown  that  C  can  be  recursively  estimated  and  used  to  backsolve  equation 
(6)  with  computational  complexity  0(M2),  This  method  converges  much  faster 
than  the  gradient  descent  method  and  is  ideally  suited  for  implementation  on 
systolic  array  processors. 

Method  5:  Dominant  Subsoace  Inversion 

If  the  application  will  allow  setting  an  upper  limit  on  the  number  of  strong 
interferences  that  the  MVDR  beamformer  is  to  reject,  then  computational 
complexity  can  be  reduced  beyond  that  of  methods  3  and  4.  Suppose  that 
K«M  represents  the  acceptable  upper  limit.  Abraham  and  Owsley^ 9  have 
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shown  that  w  can  be  approximated  by  an  expression  involving  the  K  largest 
eigenvalues  of  Rx  and  their  associated  eigenvectors  (which  are  said  to  span  the 
signal  subspace  of  Rx ),  and  Yang  and  Kaveh^o  have  shown  that  these 
quantities  can  be  computed  with  0(MK2).  This  method  is  very  well  suited  for 
implementation  on  either  systolic  or  SIMD  parallel  processors. 

Because  method  1  has  computational  complexity  0(M3)  and  method  3  is 
not  well  suited  for  running  on  parallel  processors,  only  methods  2,  4,  and  5  are 
serious  contenders  for  implementation.  The  requirements  of  any  specific 
application  should  then  dictate  which  of  the  three  methods  to  use. 


SUMMARY 


This  report  provides  an  introduction  to  the  subject  of  beamforming  for 
sonar  applications.  It  explains  why  we  need  directional  antennas,  what  the 
difference  is  between  continuous  apertures  and  arrays,  and  how  we 
characterize  performance  by  means  of  beam  patterns,  directivity  index,  and 
array  gain.  Single  beam  and  multibeam  systems  are  described  as  they  have 
evolved  from  the  1940's  to  the  present.  Also,  the  important  topic  of  adaptive 
beamforming  is  described,  including  discussions  of  various  criteria  of  optimality, 
a  null-steering  interpretation  of  ABF  performance,  a  comparison  of  beam  space 
versus  element  space  ABF,  and  five  methods  for  ABF  implementation. 

This  report  will  be  followed  with  a  classified  addendum  that  will  review 
past  efforts  to  evaluate  ABF  performance  on  data  from  actual  Fleet  arrays.  Also 
included  will  be  processing  requirements  for  real-time  ABF  implementation  with 
a  representative  towed  array. 
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